Giant spin-polarized current in a Dirac fermion system at cyclotron resonance 
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We report on the observation of the giant spin-polarized photocurrent in HgTe/HgCdTe quantum 
well (QW) of critical thickness at which a Dirac spectrum emerges. Exciting QW of 6.6 nm width 
by terahertz (THz) radiation and sweeping magnetic field we detected a resonant photocurrent. 
Remarkably, the position of the resonance can be tuned from negative (-0.4 T) to positive (up 
to 1.2 T) magnetic fields by means of optical gating. The photocurent data, accompanied by 
measurements of radiation transmission as well as Shubnikov-de Haas and quantum Hall effects, 
give an evidence that the enhancement of the photocurrent is caused by cyclotron resonance in a 
Dirac fermion system. The developed theory shows that the current is spin polarized and originates 
from the spin dependent scattering of charge carriers heated by the radiation. 

PACS numbers: 73.21. Fg, 72.25.Fe, 78.67.De, 73.63. Hs 
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The electron dc transport in semiconductor systems 
with massless Dirac fermions has recently moved into 
the focus of modern research yielding challenging fun- 
damental concepts as well as holding a great potential 
for applications [H-Q. The linear energy spectrum al- 
lows the observation of quantum kinetic effects and, on 
the other hand, gives a rise to a new class of phenom- 
ena absent in materials with parabolic dispersion. The 
massless Dirac fermions are realized in graphene [4| , at 
surface states of bulk topological insulators (TI) [JH, Q , 
in edge channels of two-dimensional TI Q as well as in 
HgTe/HgCdTe QWs of critical thickness 0-0 • In the 
latter case and TIs, the linear energy spectrum is formed 
by strong spin-orbit interaction which locks the orbital 
motion of carriers with their spins. The interest in Dirac 
fermions in such materials resulted in theoretical con- 
sideration and observation of such fundamcntalphysical 
phenomena as the quantum spin Hall effect 

Sana- 

12|, quantum Hall effect (QHE) on topological surface 



states [13j], magneto-electric effect [lj, |15|, and quan- 
tum interference effects fl6rl2ll ]. A considerable atten- 
tion has been also given to the nonlinear high frequency 
(HF) transport phenomena. A plethora of such effects 
has been treated theoretically, including photogalvanics 
in TI systems 22h25|| . and second harmonic generation 
(SHG) [H as well as radiation-induced QHE and 
topological states [H] • While a great number of propos- 
als have been published in the last two years, the num- 
ber of experiments on the topic is limited so far by a 
few publications reporting the observation of SHG and 
photogalvanic effects in 3D TIs induced by near infrared 
radiation [H, [H-[H | . 



Here, we report on the observation of a dc current ex- 
cited by THz radiation in HgTe/HgCdTe QWs of critical 
thickness. We show that the current is giantly enhanced 
at cyclotron resonance (CR) being a few orders of mag- 
nitude higher than THz radiation excited photocurrents 



detected in other non-magnetic QW structures. Due to 
the non-equidistant energy separation of Landau levels 
in systems with linear electron dispersion, the CR posi- 
tion is tuned by the variation of carrier density applying 
optical gate. The microscopic origin of the current is dis- 
cussed in terms of the cyclotron motion, spin-dependent 
scattering and Zeeman splitting. We show that the cur- 
rent is spin-polarized and its enhancement comes from 
constructively contributing three factors: strong spin- 
orbit coupling, large ^-factor in HgTe/HgCdTe QWs, and 
efficient radiation absorption at CR. 

The experiments are carried out on (013)-oricntcd 
HgTe/Hgo. 3 Cdo. 7 Te QWs j^. Single QW samples with 
widths, L w , of 6.6 nm, and 21 nm and mobilities about 
10 5 cm 2 /(V-s) at T = 4.2 K are investigated. The struc- 
tures cross section is shown in Fig.QJa). Eight ohmic con- 
tacts have been prepared at the corners and in the middle 
of the edges of 5 x 5 mm 2 samples. Magneto-transport 
measurements show well pronounced Shubnikov-dc Haas 
oscillations and QHE plateaus, see Fig. [TJc). To achieve 
a controllable variation of the carrier density we applied 
optical gating using the persistent photoconductivity ef- 
fect well known for HgTe/HgCdTe QWs jSlsllii- We 
illuminate the sample by red light emitting diode for a 
time ti resulting in a change of the carrier density (type) 
which could be restored by heating the sample above 
T ps 150 K. The carrier densities measured for different 
U used in experiments are given in Table HI 

For photocurrent excitation we apply a cw CH3OH 
laser emitting a radiation with frequency / = 2.54 THz 
(wavelength A = 118 /jm) [35| . The radiation with power 
PwlOmW is focused in a spot of about 1.5 mm di- 
ameter and modulated at 800 Hz. Right (ct + ) and left 
(<r _ ) handed circularly polarized light is obtained by a 
A/4-plate. The experimental geometry is sketched in 
Fig. H{b). In (013)-oriented QWs, excitation by nor- 
mally incident THz radiation results in a photogalvanic 
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FIG. 1: (a) Cross section of the investigated structures, (b) 
Experimental set up used for the photocurrent and transmis- 
sion measurements, (c) Magneto-transport data obtained in 
the van der Pauw geometry for 6.6 nm HgTe/HgCdTe QW 
sample without (top) and with (bottom) optical gating. The 
change of slope in the Hall signal indicates that initially p- 
type sample becomes n-type due to optical gating. 



current even at B z = 3(1, 3?| , see Supplementary Mate- 
rials. Owing low symmetry of QW, the photocurrent has 

no predefined direction and its magnitude, j 



' 3x J y 5 

can be deduced by measuring the signals along two or- 
thogonal directions [38j . The current-induced photovolt- 
ages U XyV are picked up across a 1 Mf2 load resistor ap- 
plying lock-in technique. The magnetic field B z up to 
4T is applied normal to the QW plane. The photocur- 
rent studies are accompanied by optical transmission, see 
Fig. [ljb), and magneto-transport experiments. 

We start with the data obtained on the 6.6 nm QW, 
which should have a close to linear dispersion |{J I2H ]. 
Exciting the sample with right-handed circularly polar- 
ized radiation and sweeping magnetic field we observed 
a strong resonant photocurrent at B c = —0.42 T, see 
Fig. [2] The signal at the resonance is more than two or- 
ders of magnitude higher than that detected at B z = 0. 
By changing the carrier type from hole to electron one 
the resonance jumps from negative to positive B z and 



TABLE I: Sample parameters measured at T = 4.2 K. 
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FIG. 2: Signal normalized by the radiation power U/P ex- 
cited by <7 + radiation vs. B z . The data are shown for mag- 
netic fields lower than the full quantization limit [39|. The 
inset shows data for <r + and <j~ light. 



moves towards higher field, being now for the electron 
density ri\ sa 2 x p 1 at B c = +0.69 T. Remarkably, at fur- 
ther increase in the electron density, the resonance posi- 
tion drifts to even higher B c , being 1.2 T for 712 ~ 3 x m, 
see Fig. [21 Switching the radiation hclicity from cr + to o~ 
changes the current sign and mirrors the results with re- 
spect to the magnetic field polarity, as shown for p-type 
conductivity in the inset in Fig. [5J The above behav- 
ior is observed for the temperature range from 4.2 up to 
150 K. The resonances are detected in the transmission 
measurements as well. The data for different carrier den- 
sities (type) are shown in Fig. [3] demonstrating a good 
correlation between the positions of the dip in the trans- 
missivity and the resonant photocurrent. 

Resonant photocurrents are also detected for the 21 nm 
QWs, a structure characterized by a nearly parabolic dis- 
persion. Photocurrent (see Fig. 2]) and transmission mea- 
surements (see Fig. [3]) clearly show that the resonance po- 
sition in these QWs is shifted to much higher magnetic 
fields B c « 3 T. Furthermore, the resonance field now only 
slightly depends on the carrier density 

The observed coincidence of peak positions in the pho- 
tocurrent and transmissivity unambiguously proves that 
the resonant current is caused by CR This is also 
supported by the fact that for a fixed radiation helicity 
the resonances in the photocurrent and transmissivity 
are detected for one polarity of magnetic field only. The 
striking fact is that, depending on optical gating, the 
resonance for QWs with L w = 6.6 nm is detected for neg- 
ative as well as for positive magnetic field and its peak 
position drastically depends on the Fermi level. These re- 
sults are in agreement with Ref. Q which concludes that 
the energy dispersion in similar structures is close to lin- 
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FIG. 3: Transmissivity of circularly polarized light as a func- 
tion of the magnetic field modulus, \B Z \. Full lines are fits by a 
superposition of Lorentz functions. The inset shows the data 
obtained for a + and a~ light for the L w =6.6nm samples. 



car. Such an electron spectrum lacks the band gap and, 
therefore, allows an easy transition from n- to p-type con- 
ductivity, as proved in transport and CR experiments, 
see Figs. [l]|3l Furthermore, recent study of gated Hall 
bar 6.6 nm QWs samples prepared from the same batch, 
as samples studied here, manifests a weak localization 
effects even in the vicinity of Dirac point giving an ad- 
dition al p roof for existence of gapless dispersion in these 
QWs [2l|. As a matter of fact, the cyclotron frequency 
in a system with linear dispersion depends strongly on 
the Fermi energy - a characteristic behavior observed in 
our 6.6 nm sample. Indeed, the cyclotron frequency is 
described by the well known expression ui c = \eB z \/m c c, 
where e is the carrier charge, c is the speed of light, and 
m c is the effective cyclotron mass at the Fermi energy. 
The latter, given by m c = pp/ (cLEf / dpp) with p? being 
the Fermi momentum, yields m c = E-p/v 2 for a system 
with linear dispersion characterized by a constant veloc- 
ity v. Taking Ep = \j2Ttn(Tiv) into account, we obtain 
for the CR position \B C \ = y/2nn(cfioj) /\ev\. From the 
resonance positions measured for electron densities n\ 
and 7T-2 (Fig. [2]) we find that the electron Fermi velocity 
is almost constant, being equal to 7.2 x 10 5 m/s. The 
value is in a good agreement with the electron velocity 
for 2D Dirac fermions in HgTc/HgCdTc QWs of critical 
thickness, v = 6.3 x 10 5 m/s,obtained from the energy 
spectrum calculated in Ref. [8(. The hole velocity, de- 
duced from our data for the density pi, is also close to 
this value (7.5 x 10 5 m/s). Note that the obtained values 
are close to the carrier velocity in graphene (10 6 m/s [H). 
The substantially higher resonance field in 21 nm QWs as 
well as the observed weak dependence of its position on 
the electron density correspond to the CR behavior in 
HgTe/HgCdTc QWs characterized by a nearly parabolic 
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FIG. 4: Left panel: Magnetic field dependence of the signal 
U/P excited by a + light in the 21 nm sample. Right panel 
shows the magneto-transport data. 



dispersion with large Zeeman splitting 

Now we turn to the microscopic origin of the current 
generation at cyclotron resonance. The fact that the 
current is excited at normal incidence of radiation im- 
plies that it is inherently related to asymmetry of car- 
rier relaxation (excitation) in fc-space in (013)-oriented 
QWs of low spatial symmetry [45j. In HgTc/HgCdTc 
QWs with known strong spin-orbit coupling and en- 
hanced magnetic properties, the most likely candidate 
responsible for the microscopic origin of the asymme- 
try is spin-orbit coupling, in particular, spin-dependent 
scattering. The latter is a known fact for other III-V 
heterostructures 4(| 47 1. We consider the resonant ab- 
sorption of THz radiation, which leads to the strong elec- 
tron (hole) gas heating. The resulting steady-state non- 
equilibrium carrier distribution is formed by the energy 
gain due to the radiation absorption, electron-electron 
collisions thcrmalizing the electron gas, and the energy 
loss due to emission of phonons. The matrix element 
of electron scattering by phonon contains asymmetric 
spin-dependent terms (odd in the electron wave vector), 
which are similar to the Rashba and Dresselhaus spin- 
orbit terms in the energy dispersion 01- Due to the 
spin-dependent part of the electron-phonon interaction, 
the energy relaxation of carriers in the spin cones is asym- 
metric and the relaxation rates for positive and negative 
wave vectors, say in x-direction, arc different |49| . The 
asymmetry causes imbalance in the carrier distribution 
in fc-space and, hence, electron fluxes, see Fig. EJb). The 
latter have opposite directions in the spin-up and spin- 
down cones. As, besides the cyclotron motion, the mag- 
netic field splits the spin cones due to the Zeeman effect, 
one of the cones is preferentially populated compared to 
the other. Consequently, the fluxes in the spin cones 
are unbalanced, and a net electric current emerges. This 
process is sketched in Fig. [5] and, besides the linear dis- 
persion, is alike that known for other non-magnetic QW 
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spin-down 



FIG. 5: Microscopic model, (a) Energy dispersion with spin- 
up (left, blue) and spin-down (right, red) cones shifted due 
Zeeman effect to higher and lower energies, respectively, (b) 
Non-equilibrium carrier distribution at Fermi energy caused 
by the spin-dependent scattering. Full and open circles sketch 
the electron distribution in fe-space. (c) Energy relaxation 
of an electron gas heated as a result of the CR absorption. 
Due to the spin dependent scattering, the relaxation rates 
for electrons with positive and negative fc are different. Bent 
arrows show schematically a predominant energy relaxation in 
the spin-up (left) ans spin-down (right) cones. The scattering 
asymmetry within each cone results in the oppositely directed 
fluxes j+ and j- shown by horizontal arrows. As due to the 
Zeeman effect the spin-down cone is larger populated, the flux 
j- is stronger than j+ and a dc electric current emerges. 



structures [46|, |47[. Obviously, the current magnitude 
is proportional to the radiation absorption, the strength 
of spin-orbit coupling, and the Zeeman splitting. The 
absorption is strongly enhanced in HgTc/HgCdTc QWs 
at the CR condition, resulting in the giant photocurrent 
observed in the experiments. While the models of the 
current formation for QWs with linear and parabolic dis- 
persions are similar, the position of CR and its behav- 
ior upon variation of the Fermi energy are different. In 
particular, for a fixed radiation frequency, the cyclotron 
frequency is almost independent of the carrier density in 
QWs with nearly parabolic spectrum. By contrast, for 
QWs with linear spectrum, it drastically depends on the 
carrier density and for circularly polarized light may even 
change its sign in the same structure. 

Following the model above, we now develop the quasi- 
classical theory of the effect. In this approach, the elec- 
tron fluxes in the cones are given by j s = e^ fc vSf s k, 
where 8f a k is the anisotropic part of the distribu- 
tion function to be found from the kinetic equation, 
s is the index enumerating the Dirac cones. Here, 
s ="+" corresponds to the cone formed from the 
states \El,+l/2) = fi(z)\T 6 , +1/2) + f 4 (z)\T 8 , +1/2) 
and \Hl,+3/2) = f 3 (z)\T s , +3/2), s = "-" corre- 
sponds to the cone formed from the states \E1, —1/2) = 
/!(z)|r 6 ,-l/2) + / 4 (z)|r 8 ,-l/2) and \Hl,-3/2) = 
/3(z)|r 8 ,-3/2), with /i(z), f 3 (z), and f 4 (z) being the 



envelope functions 0. 

We consider that at low temperatures, relevant to the 
experimental conditions, the momentum relaxation of 
carriers is limited by clastic scattering from static de- 
fects while the energy relaxation is governed by deforma- 
tion interaction with bulk acoustic phonons. The spin- 
dependent asymmetry of electron-phonon interaction in 
(013)-grown QWs is caused by the strain-induced cou- 



pling between the T 6 and T 8 band states [50h52| and 
deviation of the QW plane from the (001) plane by the 
angle 9 18.4°. Then, for the linear electron disper- 
sion in QWs with symmetric confinement potential, the 
resulting fluxes have the form (see Supplementary Mate- 
rial) 



J±,x = +" 



\e\v sin 29 ni 



d 



2\/2 



3±,: 



+ - 



evsm29 n± 



2\/2 



dEp 1 + (w c 7"p) 2 



dEF 1 + {u) c t p )' 



Vv, (i) 



Here a: || [100] and y||[031] are the in-plane axes, n± are 
the carrier densities in the cones, n = n + + n_, t p is 
the momentum relaxation time, / is the radiation inten- 
sity, r\ is the free-carrier absorbance, £ is a dimensionless 
parameter, 

JLoo ^cv Z\ 3 (^ C Z h + ^ V Z 33 + ^ v Zn)q 2 z dq z 
/+^[(S C Z 11 +S U Z 44 )2 + (S„Z 3 3) 2 ]^ 2 ' 

Zij = fi(z)fj(z)exp(iq z z)dz, E c and E„ are the 

deformation-potential constants in the bands Tq and 
T 8 , respectively, and 5™ is the interband deformation- 
potential constant [HJ, [52[ • 

For the circularly polarized radiation i] in the vicinity 
of CR has the form 



V 



2e 2 E F 



ctfriu 1 + (w - w c ) 2 t2 ' 



(2) 



where n u is the refractive index. Assuming that the mag- 
netic field splits the electron states due to the Zeeman 
effect but does not affect the scattering and taking into 
account that in experiment lu c t p 3> 1, we estimate the 
net electric current j = j + + j_ as 



. _ csin2gg/xo 
3 ~ 2V2v E F 71 ' 



(3) 



where (io is the Bohr magneton. Following the spectral 
behavior of the absorbance 77 the current j exhibits a 
sharp peak at CR, whose position depends on the Fermi 
energy, in accordance with the experiment, see Fig. [2] 

To summarize, we demonstrate that CR absorption by 
Dirac fermions in HgTe/HgCdTe QWs of critical thick- 
ness results in a resonant spin polarized electric cur- 
rent. The effect is very general and can be observed 
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in other Dirac fermion systems with a strong spin or- 
bit coupling, e.g., surface states in 3D topological insu- 
lators like Bi 2 Se 3 and Bi 2 Te3. For the latter case such 
a study is of especial interest, because these crystals are 
centrosymmetric and, due to symmetry reasons, the pho- 
tocurrent emerges only at surface. Consequently, it pro- 
vides a unique selective access to fine details of their band 
structure like, e.g., effective mass and group velocity, as 
well as to the spin transport and spin-dependent scatter- 
ing anisotropy [53L [H[ • Finally, large resonant currents 
detected at low magnetic fields, about 0.5 T for 2.5 THz, 
indicate that HgTc/HgCdTc QWs of critical thickness 
are a good candidate for frequency selective CR-assistcd 
detectors similar to that based on photoconductivity in 
bulk InSb (HI, but operating at about 10 times lower 
magnetic fields. 
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SUPPLEMENTARY MATERIAL 

PHOTOGALVANIC EFFECTS AT ZERO 
MAGNETIC FIELD 

Illuminating our (013)-oriented QW structures with 
normally incident THz radiation we observed a dc electric 
current even in the absence of magnetic field. Figure [5] 
shows the dependence of the photosignal U y /P oc j y /P on 
the azimuth angle a measured for 6.6 and 21 nm QWs, 
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respectively. Here a is the angle between the light po- 
larization plane and x-dircction. In both cases the data 
are well fitted by U y (a)/P = A + B sin(2a) + C cos(2a). 
While in 6.6 nm QWs such a current has not been de- 
tected so far, in 21 nm QWs it has been studied in a very 
details in Ref. [l[ and demonstrated to be due to the 
photogalvanic effect 0, 0]- A particular feature of the 
photogalvanic effect in QWs of (013)-orientation orienta- 
tion is that, in contrast to (OOl)-grown QWs, it can be 
excited at normal incidence of radiation. Quantum wells 
grown on the (013)-oriented substrate belong to the triv- 
ial point group Ci lacking any symmetry operation ex- 
cept the identity. Hence, symmetry does not impose any 
restriction on the relation between radiation electric field 
and photocurrent components. The polarization depen- 
dence of the photocurrent in structures of the Ci point- 
group symmetry for the excitation along the QW normal 
with linearly polarized light is given by [l[ 



Xxxy sin 2a 



Xxxx Xxyy 



COS 2a 



Jy 



Xyxy sin2a 



Xyxx + Xyyy 



(4) 



(5) 



Xyxx Xyyy 



cos 2a 



Ir,. 



where x is the third rank photogalvanic tensor. Ex- 
actly this polarization dependence is observed in exper- 
iment, see Fig. [51 For Ci-symmetry group all compo- 
nents of the tensor % are linearly independent and may 
be nonzero. Consequently, even for a fixed light polariza- 
tion, the CPGE photocurrent direction is not forced to 
a certain crystallographic axis. Moreover, it varies with 
temperature, radiation wavelength, etc [l[ . 
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FIG. 6: Polarization dependence of the photovoltage normal- 
ized by radiation power, U y /P, excited with normal incident 
radiation in (a) 6.6 nm and (b) 21 nm QWs at zero magnetic 
field. The normalized photosignal is shown for different car- 
rier densities and fitted after U y (a) / P = A + B sin(2a) + 
Ccos(2a). 



MICROSCOPIC THEORY 

We describe the electron states in HgTe/CdHgTe quan- 
tum wells in the isotropic k-p model following Ref. 
The Dirac cones in the QW of critical width are formed 
from the four states 

\El, +1/2) = / 1 (z)|r 6 ,+l/2> + /4(z)|r 8 ,+l/2>, 

|m,+3/2) = / 3 (z)|r 8 ,+3/2), 

|£l,-l/2) = hiz^Te, -1/2) + U(z)\T 8 , -1/2) , 

\Hl,-3/2) = / 3 (z)|r 8 ,-3/2), (6) 

which are degenerate at k = 0, with k being the in- 
plane wave vector. Here fi(z), fz{z), and fi{z) are the 
envelope functions, which can be chosen real, z is the 
growth direction, |r 6 ,±l/2), |r 8 ,±l/2), and |r 8 ,±3/2) 
are the basis functions of the T§ and Tg bands. At k ^ 



0, the states (JS]) are coupled that is described by the 
effective Hamiltonian 

/ iAk + \ 
H= -iAk- 

-iAk- 
\ iAk+ / 

where A is the (real) constant describing the in- 
plane velocity, A w (P/V§) / fi(z)f 3 (z)dz, P = 
i(h/mo)(S\p z \Z) is the Kane matrix element, and k± = 
k x ± iky. The sign of A depends on the sign of fz{z), we 
take A > 0. Solution of the Schrodinger equation with 
the Hamiltonian (fT]) for positive energy Ek = Ak yields 
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two functions 



i'+.k 



exp(ifc ■ p) 
V2 



1 

—i exp(-itp) 





(8) 



where S, X\ Y', and Z' are the Bloch amplitude of the 
Tq and Tg bands, respectively, referred to the cubic axes 
x'||[100], [010] , and z'||[001], 6 w 18.4° is the angle 
between the [001] and [013] axes, and the symbols j" and 
I denote the spin projections +1/2 and —1/2 onto the z 
axis, respectively. 



i>-,k 



exp(ifc • p) 
V2 



\ i exp(itp) / 



where ip = a,Tcta,n(k y / k x ) is the polar angle of the wave 
vector. 

We consider (013)-grown QWs and chose the coordi- 
nate frame [100] , y||[031], and z||[013]. In such coordi- 
nate system, the basis functions of the and T% bands 
can be presented in the form 



|r 6 , 
|r 6 , 



-1/2) 
-1/2) 



S J, , 



(9) 



|r 8 , +3/2> = 
|r 8 ,+i/2) = 

|r 8j -i/2> = 
+ 

|r 8j -3/2) = 



X' + i(Y' cos 9 


- Z' sin 


V2 




^(Z' cos0 + Y" 


sin 9) t 


X' + i(Y' cos0- 


Z' sin6») 






\J^{Z' cos + Y' 


sin 6) I 


X'-i(Y' cos 9- 


Z' sin0) 






X'-i(Y' cos 9- 


Z' sin 9) 


V2 



t, 



t, 

I, (10) 

I 



The deformation interaction of electrons with acous- 
tic phonons in zinc-blende-type crystals has both intra- 
band and interband contributions [H, @] . The matrix ele- 
ments of strain-induced interband coupling are given by 

Vs,X' = ^cvUy'z', Vs y Y' = ScUx'y" Vs,Z' = EcvUx'y', 

where S cl) is the interband constant of the deformation 
potential and u a fs are the strain-tensor components used 
here in the primed coordinate system. Note, that E cv / 
in non-centrosymmetric crystals only. The matrix ele- 
ments of strain-induced intraband interaction are taken 
in the form Vs,s = S c Trit a/3 , V x >,x< = Vy>,y> = V z >,z> = 
E„ Trw a ^, where S c and H„ are the deformation-potential 
constants (3„ = a in the Bir-Pikus notation Q, the con- 
stants b and d are neglected for simplicity) . Accordingly, 
the Hamiltonian of electron-phonon interaction in the ba- 
sis of functions © and (ITU1) has the form 



V 




(11) 



where V c = E c (Tru a /3)h, V v = S„(Tr u a p)h, I 2 and I 4 
are the identity matrices 2x2 and 4x4, respectively, 



/ (u yz — iu xz ) cos 29 + (u zz /2 — u yy /2 + iu xy ) sin 29 

71 



cv ~ — Ct) 



(u xy cos 29 + u xz sin 29) 



{uyz + iu xz ) cos 20 + (u zz /2 — Uyy/2 — iu xy ) sin 29 

o 





(M yz - iu xz ) cos 26> + (u zz /2 - u yy /2 + iu zl) ) sin 29 

— (u^y cos 29 + u xz sin 29) 
(u yz + cos 26* + (u zz /2 - u.yy/2 - iu xy ) sin 20 
V2 



r 



and the strain-tensor components u a p are rewritten in 
the QW coordinate frame. 

The dominant contribution to electron scattering is 



given by the terms proportional to u zz because the out- 
of-planc component q z of the wave vector of the phonon 
involved is typically much larger than the in-plane com- 
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ponent gii. In this approximation, the matrix element is limited by elastic scattering from static defects while 

of electron scattering from the state (s,k) to the state the energy relaxation is governed by deformation interac- 

(s,fc'), described by the wave functions ©, assisted by tion with bulk acoustic phonons. Accordingly, 5f s k can 

emission or absorption of a bulk acoustic phonon with be found from the Boltzmann equation 
the wave vector q has the form 

12 .,. , ^[«*fl]-^ = Er^&/"( 1 -/. fc ) 



y(±) 
sk' ,sk 



-I 



<±) 



hN tM 

2pfl q 



i sin 28 



2\/2 



where 



'^cv (c 



No 



isip 



-isip 



)Zi 



(12) 



(1 ± l)/2, N q is the phonon oc- 



cupation number, p is the crystal density, D, q = ciq is 
the phonon frequency, q is the speed of sound, Z^ = 

J-™ fi( z )fj( z ) exp(iq z z)dz, and s ="±" is the index 
enumerating the Dirac cones. The matrix elements ([12")) 
contain asymmetric terms which are responsible for the 
emergence of oppositely directed electron fluxes j s in the 
cones during the energy relaxation of heated electron gas. 

To calculate the electron fluxes, we introduce the elec- 
tron distribution function f s k = / s fc + Sf s k, where f s k is 
the quasi-cquilibrium function of the Fermi-Dirac type, 
6 f s k is the anisotropic part of the distribution function. 
It is assumed that the radiation absorption followed by 
electron-electron collisions forms the quasi-equilibrium 
electron distribution with the electron temperature T e 
which is slightly higher that the crystal lattice temper- 
ature To. The electron temperature can be found from 
the energy balance equation 



E 

s.k.k' 



w. 



(Ph) 



sk' .8 



k (e k - Sk>)fsk(l - fsk') = Ir) 



(13) 



where Wj$„ k 



(VR) E,,± \v£l k \*8( Ek , - e k ± nn q ) 

is the rate of electron scattering assisted by a phonon 
emission and absorption, scattering processes between 
states in different cones are neglected, / is the radia- 
tion intensity, and ry is the free-carrier absorbancc. The 
left-hand side of Eq. (fT31 describes the electron energy 
losses due to cooling by phonons while the right-hand side 
stands for the energy gain by the free-carrier absorption 
of radiation. 

The electron fluxes are determined by the anisotropic 
part of the distribution function 



k 



(14) 



where v = Vfee/?i = vk/k is the electron velocity. We 
consider that at low temperatures, relevant to the exper- 
imental conditions, the momentum relaxation of carriers 



w^Lf sk (i-f sk ')}- 



5 fsk 



(15) 



where t p is the momentum relaxation time. The straight- 
forward calculation shows that the electron fluxes in the 
cones have the from 



3±,x = T- 



\e\v sin 28 n± I d ^cT v 



2y/2 n \ dEp 1 + {uj c T p f 



(16) 



ev sin 26 n± ( d t„ 

J± - = T ^71"~fei + (^) 2)e/?? ' (17) 



where n± are the carrier densities in the cones, n = 
n + + ri- is the total density, tu c = \e\v 2 B / (cEp) is the 
cyclotron energy, and 



S-oo ^cvZn^cZn + B„Z33 + z. v Zi4)q 2 dq 
/+~[(5 c Z n + E v Z 44 y + {~ v Z^f]qHq 



Equations (fl~6)) and (jTTj) are obtained assuming that 
Ep ksT hQ q and |n+ — n_| <C n. 
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